Robustness of the optical-conductivity sum rule in Bilayer Graphene 
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We calculate the optical sum associated with the in-plane conductivity of a graphene bilayer. 
A bilayer asymmetry gap generated in a field-effect device can split apart valence and conduction 
bands, which otherwise would meet at two K points in the Brillouin zone. In this way one can 
go from a compensated semimetal to a semiconductor with a tunable gap. However, the sum rule 
turns out to be 'protected' against the opening of this semiconducting gap, in contrast to the large 
variations observed in other systems where the gap is induced by strong correlation effects. 



I. INTRODUCTION 

Sum rules on the conductivity have played an impor- 
tant role in the analysis of optical conductivity data, 
which give information on electron dynamics. In its sim- 
plest form, for an infinite free electron band, the sum 
rule gives the plasma frequency which is independent of 
temperature, interactions and impurity scattering. 1,2 In 
tight-binding models it is related to the second derivative 
with respect to the momentum k x , v of the band disper- 
sion £k, times the probability of occupation of the state 
|k). 1,2 In the special case where only nearest-neighbor 
hopping is present on a square lattice, the sum-rule inte- 
gral reduces to minus one half the kinetic energy in ap- 
propriate units. As a consequence, in general it is affected 
by the interactions present in the system, and these de- 
termine both its absolute value and its variations with 
temperature. This fact has been investigated recently 
both experimentally and theoretically in the context of 
high-T c cuprate superconductors (see Ref. [1,2] for a re- 
view). It follows from the Sommcrfcld expansion that the 
relative temperature variation of the sum-rule integral 
in non-interacting tight-binding models is proportional 
to (T/t) 2 , where t is the nearest- neighbor hopping. In 
the experiments done in the normal state of cuprates, 
the temperature variations of the the sum rule are much 
larger than this estimate, and in some cases deviate from 
a simple T 2 law. The discrepancy can be attributed to 
correlation effects. For example, in the studies of the 
Hubbard or t — J model 3,4 , a new energy scale emerges 
associated with the reduced width of the renormalized 
band. This parameter replaces the hopping parameter 
t in the {T/t) 2 dependence leading to its enhancement. 
It is shown in other theoretical approaches that the sum 
rule integral measures instead a specific average of the 
quasiparticle lifetime 5-8 normalized to the Fermi energy 
ef, and hence such optical experiments ultimately probe 
correlation effects through lifetime broadening. These ex- 
amples illustrate that optical data can provide important 
insight into correlation effects. 

Graphene, which is a single layer of graphite, has re- 
cently been isolated and its properties investigated (see 
Ref. 9 for a review). Bilayer graphene 10 and thicker 



graphite films are also now widely produced. Since 
graphene possesses truly remarkable properties both 
from the technological and theoretical points of view, 
there has already been considerable work done on mono- 
layer and bilayer graphene and also on related materials. 
In particular, recent data arc available on the AC conduc- 
tivity in the infrared region of graphite 11,12 , several layer 
epitaxial graphite 13 , and on monolayer 1 and bilayer 15 
graphene. There has also been much theoretical work on 
the microwave and infrared conductivity of graphene 1 
(see Ref. 19 for a review) and on a bilayer 20-22 and 
multilayer. 18,23 Thus, it is of interest to anticipate some 
general results on the optical-conductivity sum rule be- 
havior of these systems, which are likely to be tested ex- 
perimentally in the near future. The optical sum rules for 
the (in-plane) longitudinal and Hall conductivities were 
studied for monolayer graphene in Ref. [24] , where the un- 
usual dependences of the sum rules on temperature and 
chemical potential were revealed. For example, at the 
neutral Dirac point the linear dispersion law of quasipar- 
ticles leads to a T 3 dependence of the sum-rule, instead 
of the T 2 law found in tight-binding bands 1,2 . 

A specific and very useful feature of graphene, not 
available in ordinary metals, is that it is possible to 
change the chemical potential [i and thus the number of 
carriers by tuning the gate voltage V g (V g oc [i 2 for mono- 
layer graphene) in a field-effect device. It was shown 24 
that in the limit T = the concentration-dependent part 
of the sum rule goes like (|^|/i) 3 . For a finite temperature 
T <C |/x|, there is a temperature dependent correction to 
the sum of order (\^\/t){T /t) 2 which is very much as in 
the more ordinary T 2 -law case discussed above. 

In this paper we generalize the previous work to the 
case of bilayer graphene. In the case considered most of- 
ten, a bilayer graphene consists of two coupled hexagonal 
lattices. The inequivalent sites Al, Bl and A2, B2 on the 
bottom and top graphene sheets are arranged according 
to Bernal (A2 — Bl) stacking: every Bl site in the bottom 
layer lies directly below an A2 site in the upper layer. The 
strongest inter-layer coupling between pairs of A2 — Bl 
orbitals, 71, changes the electron dispersion from a linear 
to a quadratic form 21,25 . In addition to this effect, bilayer 
graphene is the only known material in which the elec- 



2 



tronic band structure can be changed significantly simply 
by applying an electric field perpendicular to the layers. 9 
Indeed, while the unbiased bilayer system is a semi-metal, 
with parabolic valence and conduction band touching 
each other at the neutrality point, the system becomes a 
semiconductor when the two layers are electrostatically 
inequivalcnt. The corresponding semiconducting gap A 
between valence and conduction bands is tunable contin- 
uously from zero to r* 0.2 eV. 9 ' 26,27 As shown experimen- 
tally, the asymmetry gap A and the carrier concentration 
can be controlled independently through an applied gate 
voltage and chemical deposition of potassium or other 
atoms on the vacuum side of the structure. 26 ' 27 Recently, 
similar results have also been obtained in a double-gated 
structure made of bilayer graphene sandwiched in be- 
tween two gate electrodes 28 . This effect has been under- 
stood by means of both tight-binding calculations 26,29 ' 30 
and ab initio calculations 31 . 

Beyond the clear potential impact that this effect has 
on device applications, it is extremely interesting from 
the theoretical point of view to understand how this 
field-induced band transition will affect the conduction 
properties of the system. In light of the previous dis- 
cussion concerning the physics of cuprates, one could 
ask if the opening of the gap due to the metal-insulator 
transition in Hubbard-like models, or across the metal- 
superconductor transition in superconducting models, 
share commonality with the physics of bilayer graphene. 
As we shall demonstrate below, these phenomena are 
drastically different. Indeed, we show that the sum rule is 
essentially protected against the opening of the semicon- 
ductor gap in this bilayer system. Despite the fact that 
A can be as big as a significant fraction of the Fermi 
energy, the relative sum-rule changes from the unbiased 
to the biased system are of order of 10 -3 . Thus, even 
though a general redistribution of spectral weight is ex- 
pected due to the opening of the gap, the total sum rule 
is practically constant. One would have expected the 
large band structure changes involved to lead to corre- 
sponding significant changes in the optical sum. Here, 
however, we show that this is not the case. The relative 
changes in optical spectral weight induced by the opening 
of the gap must largely compensate each other, leaving 
the total weight almost unchanged. 



II. DERIVATION OF THE SUM RULE 

For a generic electronic model described by the Hamil- 
tonian H the optical-conductivity sum rule is given by: 1,2 



of H evaluated at finite vector potential A: 
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(1) 



where V is the unit-cell volume and t xx is the so-called 
mass tensor, which appears in the second-order expansion 



ff(A)«H(0)- 



-A a (n)j a (n) 



2c? 



A 2 a (n)T aa {n) 



(2) 

and j a is the current density in the a direction. The vec- 
tor potential A is introduced in the tight-binding Hamil- 
tonian (3) written below by means of the Peierls sub- 
stitution, 4n fe 2n+<S' -> 4n CX P (~ it Jn+6' Arfl ") 6 2n+5' • 

Here a2n and &2n+<5' are the Fermi operators of electrons 
on A2, B2 sublattices of the second layer. (Note that 
this Peierls substitution corresponds to the 2nd layer, for 
the 1st layer it is similar. The spin label is omitted.) 
The positions of A2 and Bl atoms are denoted as n and 
they are connected to their nearest neighbors on B2 (Al) 
sites by the three vectors 5' (5). Accordingly, the general 
tight-binding Hamiltonian for a bilayer reads 



H 



= -^( a in+<S 5ln + k - C -) ~ *X^ a 2nW«' + h.C.) 
11,(5 n,5' 

+ 71 X)( a 2n & ln + h.C.) + 73 ^2(af n b 2n +S' + h.C.), 
n n,5' 

- iA^(a+ +i a ln+<J + 6+6 ln -(1^2)), (3) 



where t is the in-plane hopping parameter, and 71 the 
strongest inter-layer coupling. When A = this Hamil- 
tonian is equivalent to the Slonzewski-Weiss-McClurc 
(SWM) model 32 ' 33 for the bulk graphite, provided that 
one retains only in-plane and 71 out-of-plane hopping, 
and assumes no dispersion along the c axis (perpendicu- 
lar to the planes) 34 . As mentioned above, the connected 
sites A2 and Bl lie directly below and above each other, 
so that the magnetic field perpendicular to the bilayer 
does not affect this coupling. This is not the case for the 
weaker Al — B2 coupling 73 included in the Hamiltonian 
(3) for completeness. In what follows, however, we will 
neglect this term, whose role is negligible as far as the 
longitudinal optical sum rule is concerned. Finally, the 
parameter A represents the asymmetry of the on-site en- 
ergies on the two layers, and it is induced by the electric 
field applied perpendicular to the graphene bilayer. As 
we shall see below, it is responsible of the splitting of the 
valence and conduction bands. 

The Hamiltonian (3) describes four bands 21 ' 25 ' 26 ' 30 
«4 a) (k), a = 1,2 with 



-(") 



(k) 



r = 



± 



|<Hk)| 2 ( 7l 2 + A 2 ) 



(4) 



where the function </>(k) = — 



ik<5, 



,-ik<5' 



By expanding 0(k) around the two inequivalent K,K' 
points in the Brillouin zone one gets the usual linear 
dependence, |</>(k)| w hvp\k\, where the wave- vector k 
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is measured from K points and vp = \/3ta/(2h) is 
the Fermi velocity, with a the lattice constant (V = 
\/3a 2 /2). Using this form, one can easily see from Eq. (4) 
that at large momenta the four bands reproduce the 



two Dirac cones \e 



Mi 



HvFk = \p\ of each uncou- 
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FIG. 1: (Color online) Top: bilayer band dispersion for A = 
(a) and A = 0.4 eV (b) as a function of p = hvpk. Bottom: 
general scheme of the experimental setup for a tunable-gap 
bilayer device (notation defined in the text). 



function (GF) G = -(T T ** t ), where we introduced 
= (atk' b 2k' a ik^ b tk)- Tnen tne averages (a+ k 6i k ) 
and (a^" k 62k) are 



pled layer. The low-energy band dispersion is, how- 
ever, drastically affected by the interplane hopping 71 
and the asymmetry parameter A. Indeed, at A = 
for \s±\ < 7i/4 the two low-energy bands are parabolic 

with the dispersion e± « ±h 2 k 2 /2m, where the effec- 
tive mass m = ji/(2v F ), and touch each other as in 
a semimetal (see Fig. la, where we used t = 3.1 eV, 
71 = 0.4 eV 26,30 ). However, when A ^ the system be- 
comes a semiconductor with a gap A = A71 / \J ^\ + A 2 „ 
at the momentum p m /hvF, that corresponds to the en- @ (i w n)k) 
crgy p 2 m = (A 2 /2)( 7l 2 + A 2 /2)/( 7l 2 + A 2 ), see Fig. lb. 



with z = iu 



G41 (w n ,k), 



(6) 



where iuo n is the fermionic Matsubara frequency. The 
corresponding elements G41 and G23 of the GF G can be 
found from the inverse GF 21 



/2 + iA -0*(k)\ 

z-|A -cj>(k) 
-</>*(k) z-§A -71 
V -0(k) -71 z + ^Aj 



/x. Then 



(7) 



r k\ #)H#)P + (A/2^) 2 ] m 

0411 [z2_ (£ (l))2][ z 2_ (£ (2) ) 2] « W 

and G23 = G4i(— A). Calculating the Matsubara sum in 
Eq. (6) and using the identity cf>(k)(d 2 (j)* (k) / dk^) + c.c = 
-(a 2 /3)|(/>(k)| 2 , we finally obtain 



2a 2 



V 



3h 2 VN 



where 



Y, i0(k)i 2 /(4 Q) (k)- M )M(4 Q) (k)), 

(9) 



k,A=±,a 



M(4">(k)) 



2 v ; 4T 



.(«) 



(k) 



(10) 



and /(e) = l/[exp(e/T) + 1] is the Fermi distribution. 
One can verify that for 7! = A = Eq. (9) reduces to 
the doubled one-layer sum of Ref. [24] . 

The expression (9) does not show any clear formal re- 
semblance to the kinetic-energy density of the bilayer sys- 
tem, that one would define using the bands (4) as: 

K = VN £ 4 Q) (k)/(4 Q) (k)-M). (11) 



k,A=±.a 



Although the linear approximation for 4>(k) is suffi- 
cient for the analysis of the low-energy band structure 
and for the numerical study of the mass tensor (9) be- 
low, the derivation of this tensor must be done using an 
uncxpanded expression for </>(k). 24 It follows from the 
definition (2) that (the factor 2 accounts for the spin) 



h 2 N 



E 



(aik&ik + a+k&ak)- 



'(k) 



dkl 



. (5) 



where we kept the full momentum dependence of 0(k) 
(N is the number of unit cells). The thermal average 
(i~xx) is calculated from the the imaginary time Green's 



Nonetheless, as we shall see these two quantities show 
approximately the same doping dependence, satisfying 
again the general relation between sum rule and kinetic 
energy, despite the complicated band evolution with dop- 
ing of the bilayer system. 



III. SUM RULE AS A FUNCTION OF DOPING 

As we already mentioned, it was shown 
experimentally 26-28 that one can tune independently 
the asymmetry gap A between the two layers and the 
total carrier density (the value of fi). This corresponds 
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to controlling the excess carrier density ni. 2 (difference 
between the densities of electrons and holes) in each 
layer as schematically shown in Fig. lc. The charge on 
the surfaces below the bottom layer and above the top 
layer is n g and no, respectively. In the experimental 
configuration of Ref. 26 n g is varied by changing the 
gate voltage and no by chemical deposition of dopants 
on top of the upper layer (see also Ref. 27), while in the 
device of Ref. 28 both no and n g are varied, by using two 
independent gates. The resulting asymmetry A between 
on-site energies in the two layers can be determined by 
equating the voltage difference A/e between the plates 
of the capacitor to its value evaluated from the charge 
density e(n + n 2 — n g - ?ii): 21 ' 26 ' 30 



A= 2T [no 



n 2 (A) 



m(A)] = 



riged n 2 (A)ed 



so 



(12) 

Here d is the bilayer thickness, £0 is the permittivity of 
free space and in the last equality we used that the to- 
tal carrier excess n = m + n 2 = —(no + n g ). Since 
m(A = 0) - n 2 (A = 0) = 0, we obtain that A = at 
n = —2tiq. The gap closes when the system is doped 
away from half-filling, as indicated by the ARPES mea- 
surements of Ref. 27. However, if the value of n is tuned 
to be exactly at zero 28 , the gap closes exactly at the neu- 
trality point, where the semimetal band structure is re- 
produced. For finite A the excess carrier densities ni i2 
are given by 

n x (A, M ) = ^J2 '(k) - M) - M Q) (k) + M) 



NV 



+fl w (k) [/(e w (k) — n) — f(e^(k) (j,)] 



(13) 



where 5 ( a )(k) = (A/£ ( K Q) (k))(l/2 + (-l) a \(j)(k)\ 2 /T), and 
n 2 (A, n) = n\{— A,/i). To determine selfconsistently 
the dependence A on n, one has to solve simultaneously 
Eq. (12) and the equation ni(A, fj,) + n 2 (A,/i) = n for 
the chemical potential fi. Once A(n) and fi(n) are de- 
termined, one can compute (t xx ) and K as a function of 
n. 

Following Ref. 27 we use d = 3.4 A, and show results 
for no = 0,77-0 = ±12 x 10 12 cm~ 2 . This choice is sug- 
gested by the fact that in the measurements of Ref. 27 the 
gap closes at n = — 2no given approximately by 24 x 10 12 
cm~ 2 (notice that n expressed in units of 10 12 cm -2 cor- 
responds to 1.9 x 10 3 n uc , where n uc is the number of 
electrons per unit cell n uc ). The doping dependence of 
the (t xx ) at T = is shown in Fig. 2. As one can see, 
even though the charging gap A changes consistently 
with doping (see inset of Fig. 2), due to the screening 
effects of the bilayer device, nonetheless the sum rule 
attains overall variations of ~ 0.1%, which are unde- 
tectable from an experimental point of view. In other 
words, despite the large optical-conductivity spectral- 
weight redistribution associated with the opening of the 



gap through the semimctal-scmiconductor transition in- 
duced by doping, the sum rule is not affected. It follows 
the same behavior that one would obtain in (electrostat- 
ically) uncoupled layers where A = at all doping (see 
the dashed-dotted line in Fig. 2). In the regime where 
the gap is relatively small (for example at negative dop- 
ing for the no = 12 x 10 12 cm -2 curve in Fig. 2) the sum 
rule follows approximately the same doping dependence 
found in the single-layer case, i.e. \T xx (n) — t xx (0)]/V ~ 
-2a 2 \^\ 3 /(9Trh 2 v 2 F ). Since 2a 2 /(97rft 2 7; 2 ) « 10- 2 cV~ 2 , 
and |/_t| 3 has an overall variation of at most O.l(eV) 3 in 
the considered doping range, one can easily get the small 
sum- rule variation with doping. At larger gap values, 
T xx (n) — t xx (0) is slightly larger than expected in the 
single-layer case, but is still too small to be detected 
experimentally As far as the temperature dependence 
is concerned we always found a (T/t) 2 variation, as ex- 
pected due to the parabolic structure of the low-energy 
bands. Including other hoppings such as 73 in the Hamil- 
tonian (3) can change the dispersion curve as discussed 
in 29 (see also Ref. 34) for energies less than 2 meV. While 
this can slightly affect the overall value of the optical in- 
tegral and of the self-consistent A, it does not change 
the main conclusion of our calculation, which focus on 
the relative doping changes of the sum rule when a large 
semiconducting gap A of order of a few hundred meV 
is opened. Having established that the relatively large 
band-structure changes caused by the opening of the gap 
A lead to a negligible changes in the optical sum, we ex- 
pect that the above mentioned hopping terms introduce 
negligible corrections as well. 
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FIG. 2: (Color online) Doping dependence of the T = sum 
rule for a bilayer system with a doping-dependent gap A(n), 
shown in the inset. We show results for different values of 
no, given in units of 10 12 cm -2 . For comparison, we show the 
equivalent doping dependence of the sum rule when no gap is 
present between the two layers (dashed-dotted line). 

Analogously, the kinetic energy, shown in Fig. 3, is not 
much affected by the opening of the gap, and it attains 
a value which is approximately 6 times the sum rule, 
t/V ~ (1/6)|AT|. There are however two differences in 
the doping dependence of the sum rule with respect to 
the kinetic energy. First, we notice that the very small 
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changes in the sum rule and kinetic energy have the op- 
posite trend: indeed, they both decrease when the gap 
opens, in contrast to what is found for example at the 
transition between a normal metal and a superconductor, 
where a kinetic-energy increase corresponds to a sum-rule 
decrease and vice-versa. This effect disappears when we 
set 71 = A = in our numerical work: indeed, in this 
case t xx reproduces exactly twice the single-layer value 
computed in Ref. [24], and the sum rule and the kinetic 
energy track each other exactly as the doping is varied. 
Second, it is clear from Fig. 3 that the chemical-potential 
jump at the neutrality point in the biased case (see in- 
set (b) of Fig. 3) is much more effective on the kinetic- 
energy doping dependence. The two curves for \K\ at 
no = ±12 x 10 12 cm~ 2 in Fig. 3 show a kink at n = 
which is clearly due to the /i(n) discontinuity at the neu- 
trality point. Indeed, in the case where no gap is present 
in the system u(n) goes smoothly across n = 0, as does 
the kinetic energy. Instead, the sum-rule curves in Fig. 
2 are not much affected by the fj, discontinuity at n = 0. 




-10 10 

n (10 12 cm" 2 ) 

FIG. 3: (Color online) (a) Doping dependence of the T — 
kinetic energy for a bilayer system with a doping-dependent 
gap A(n) and without it. Observe the small variations of K 
in the considered doping range. Insets: doping dependence of 
the chemical potential when a gap is present (b) or not (c), 
for parameters as in the main frame. 



tronic excitations, the negligible dependence of the sum 
rule on the carrier density was not a priori expected. 
Indeed, despite the profound band-structure changes in- 
duced by the change of doping, no significant signature 
appears on the overall spectral weight. This implies that 
the spectral weight lost in the gap must appear above 
it, and the two must largely compensate for each other. 
Observe that in order to be able to test experimentally 
such a prediction one would need to define an "experi- 
mental" cut-off. Indeed, even though the relation (I) is 
theoretically established by integrating the optical con- 
ductivity to all frequencies, in practice an intrinsic cut-off 
is provided by the frequency uj c above which transitions 
to other electronic bands not considered in the present 
study would develop. 2 According to optical studies on 
graphite 35 , the optical sum rule saturates to the value 
corresponding to one electron per atom (as due to the ir 
band considered here) around 8 eV, while optical tran- 
sitions coming from the a bands appear around 15 eV. 
However, in the present case a much smaller cut-off could 
be used, if one wants just to compare the spectral- weight 
variations induced by the opening of the gap. Indeed, 
the spectral-weight redistribution is expected to fall in a 
range of frequencies of order 2 A, as confirmed recently in 
Ref. 22 . Thus, an upper cut-off of order of oj c ~ 1 — 2 eV 
should be enough to account for the gap-opening effects, 
and to test experimentally the predicted robustness of 
the optical sum rule. 

It is worth noting that in our calculation only the elec- 
trostatic interaction between layers was included, since 
this is believed to be the most important effect in the 
system. However, as in the case of cuprates supercon- 
ductors other mechanisms (electron repulsion or electron- 
phonon interactions) could be at play and modify the sum 
rule. Thus, the experimental verification of our predic- 
tion would help understanding if other interactions need 
or need not be taken into account in dealing with these 
systems. 
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